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The research proposed for the third year of project was to study the structure of 
turbulence in boundary layers on aerodynamic surfaces and to develop stereoscopic PIV 
methods for use in such studies (c.f. proposal dated 9/18/95) 

This report describes the final results of the project, which are reported in the reference in 
Appendix A. This reference contains the technical details, and will constitute the main body 
of the report. The paper in Appendix A was first presented at the 8th International 
Symposium on Applications of Laser Techniques to Fluid Mechanics, July, 1996. This 
paper was one of 30 papers selected from more than 200 papers presented at the symposium 
for publication in a book by Springer -Verlag. It will become available during 1997. 

A stereoscopic PIV system was constructed using two Kodak Megaplus 4.2 CCD cameras to 
provide the maximum resolution available with CCD’s. The cameras are separated by a 30 
degree angle and the separation is purely lateral. The field of view is approximately 80mm x 
80mm at a working distance of 1 600mm. The air flow is seeded with 1 micron oil droplets 
and illuminated by synchronized YAG laser with ~ 100 mJ per pulse. The analysis of the 
CCD images is performed by software from TSI, Inc called Insight. This software was 
specially modified to provide acquisition of the data from the two CCD cameras 
simultaneously. The cameras were registered physically to within several pixels in the image 
plane, then registered electronically to correct for residual offset. Once vector displacement 
fields were acquired, they were combined using software developed in-house. The source 
code for this software is given in Appendix B. 

The flow investigation was performed using the stereo particle image velocimeter to 
provide views of the instantaneous , two-component velocity fields in streamwise-wall 
normal planes of a boundary layer in a wind tunnel. The results confirm the formation of 
near-wall shear layers, a major new insight into the structure came in the form of strong out- 
of-plane motion along the inclined layers, indicating strong rotation, even of the larger 
scales. 
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The flow investigation was performed using the stereo particle image velocimeter to 
provide views of the instantaneous , two-component velocity fields in streamwise-wall 
normal planes of a boundary layer in a wind tunnel. The results confirm the formation of 
near-wall shear layers, a major new insight into the structure came in the form of strong out- 
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Structure of a Turbulent Boundary Layer Using 
a Stereoscopic, Large Format Video-PIV 


Z. C. Liu 1 , R. J. Adrian 1 . C. D, Meinhart*- and W. Lai^ 

* Department of Theoretical and Applied Mechanics 
University of Illinois, Urbana, Illinois 61801 

^Department of Mechanical and Environmental Engineering 
University of California, Santa Barbara, California 93106 

^TSI, Inc.. St. Paul. Minnesota 55164 

Abstract. Development of a stereoscopic particle image velocimeter tor the 
measurement of three-dimensional vectors on a planar domain is described. The 
camera is based on two large format (2k x 2k) video cameras. Experiments in a 
turbulent boundary layer at Reg = 2525 demonstrate its ability to measure three- 
dimensional turbulent flow. In addition to the quantitative value of the out-ot- 
plane component, it is found that having the complete three-dimensional vector 
also significantly improves the qualitative visualization of the flow. 

Keywords, particle image velocimetry, boundary layer, three-dimensional, 
stereoscopic, videographic, turbulence 

1 Introduction 

The extension of particle image velocimetry to measurement ot three-dimensional 
velocity vector fields is both desirable and achievable. Hinsch (1995) reviews 
many of the techniques used for three-dimensional measurements. In general, there 
are two classes of measurement techniques: those which measure three-dimensional 
velocity vectors on full three-dimensional domains , i.e. volumes, and those which 
measure three-dimensional vectors on planar domains, i.e. light sheets. Three- 
dimensional volumetric measurements can be performed by using holographic 
techniques (Barnhart et al. 1994, Meng and Hussain 1994), by using 
photogrammetric particle tracking techniques (Maas, Gruen and Papantonoiu 1993, 
Brodkey 1977, Nishino and Kasagi 1991), or by scanning a laser light sheet 
rapidly ( Brucker 1995). Photogrammetric methods use three or four cameras to 
view particles from several different directions. Volumetric measurements are in 
general difficult and require equipment that is rather different than the standard 
planar PIV equipment which uses a laser light sheet and a single camera. 
Holographic systems require a unique set of techniques and apparatus, and three- 
dimensional particle tracking systems involve multiple cameras and specialized 
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software, plus they place demands on the optical access required to view the How 
field. 

However, the most demanding aspect of volumetric measurements is the 
sheer number of velocity vectors that are obtained from the measurements. A 
modest 100 x 100 x 100 velocity grid measurement yields one million vectors, 
and if these measurements are to be repeated hundreds or even thousands ot times 
for the purpose of statistical averaging, the amount of data is overwhelming. 
Thus, while volumetric measurements can be extraordinarily valuable tor studying 
instantan-eous flow field structures, they are often too rich in data to be used on a 
routine basis with present equipment. 

An intermediate approach to three-dimensional volumetric measurements 
with PIV is to perform measurements of a three-dimensional velocity vector tiled 
on a planar domain. In this way, the number of vectors is the same as in ordinary 
two-dimensionai planar PIV, but the third component contributes significantly to 
the experimenter's capability to visualize the How and it is valuable tor the 
purposes of quantitative analysis of the How. Further, tull measurement ot the 
three-dimensional vectors eliminates the perspective error that is inherent in 
monocular PIV systems (Adrian 1991, Prasad and Adrian 1993). This error can be 
quite significant if the out-of-plane velocity component has non-negiigible 
magnitude relative to the in-plane components and/or the angular held-ot-view is 
not small. 

Systems operating with this capability can be designed using various 
aspects of image correlation ( Robinson and Rockwell 1993), or using the change 
in image magnification (Willert and Gharib 1992). However, the most common 
technique is the well-known stereoscopic method whereby one obtains a pair ot 
image planes, each image plane viewing the particles in the illuminated light sheet 
plane from a different direction. Displacements of the particle images in the two 
different views differ because of the different viewing angles, and measurements of 
these displacements can be used to solve for the full three-dimensional 
displacement of particles. This technique can be applied to individual particles, as 
in particle tracking systems, or it can be used in combination with correlation 
techniques that measure the displacements of groups of particles, or even 
continuous grey level patterns. The correlation approach is desirable because the 
three-dimensional velocity vectors can be obtained on unitorm grids, and the 
problem of matching individual pairs is solved automatically. The limitation of 
the stereo technique is that to assure good measurements of the out-of-plane 
component of velocity, the viewing angles between the stereo lenses must be 
substantial, of the order of thirty degrees, and the time between exposures of the 
images must be small enough that the particles remain within the thickness of the 
light sheet. These requirements have been explored amply and shown to be 
achievable and not inconvenient in practice (Prasad and Adrian 1993, Troy and 
Adrian 1996)). 

Stereographic systems of conventional form have been used by Arroyo 
and Greated (1991), Prasad and Adrian ( 1993), and Troy and Adrian (1996). These 
experiments each used two photographic recordings, interrogated the recordings to 
obtain displacements and then solved the stereo equations to obtain three- 
dimensional velocity vectors. The difficulty with using photographic film is the 
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problem of registering the two pieces of film with respect to each other so that 
measurements of displacements locations in the first image can be identified 
precisely with locations in the second. This problem must be faced for each pair 
of photographic exposures that are taken in a stereo-photographic system. The 
registration process is more laborious and time consuming than developing the 
film in itself, and it places a fundamental constraint on the utility of this method. 

To make feasible extensive quantitative analysis of three-dimensional 
vector fields, it is necessary to make the image acquisition analysis process easy 
enough and fast enough to permit thousands of images to be taken for averaging 
purposes. To this end, the photographic recording must be eliminated, and 
replaced by videographic recording. Video cameras in a stereo system can be 
registered during the construction of the stereo camera and by this means they can 
be aligned once and for all. thereby eliminating the need for any future efforts to 
register the two images. The registration can be accomplished by means of 
mechanical alignment during construction, or by electronic image processing 
alignment during the experimental stage. The electronic alignment is necessary in 
situations where an aberrating medium is placed between the object plane and the 
cameras, so that the registration that was originally achieved in the construction of 
the camera is distorted by the medium. In this case, a calibration procedure can be 
used to correct for aberrations, provided they are not so severe as to completely 
destroy the images of the individual particles. 

In this paper, the design and construction of a stereographic camera based 
on video recording is described, and some of the procedures used in the aberration 
calibration and correction process are explained. The system is applied to the 
visualization of three-dimensional vector field structure in a low Reynolds number 
turbulent boundary layer. 

2 Stereoscopic PIV System 

The stereoscopic camera is shown schematically in Fig. 1. The camera lenses L\ 
and L 2 lie in a plane parallel to the object plane defined by the laser light sheet. 
They are offset laterally by du and d/j, respectively. Since the camera arrays also 
lie in a plane parallel to the object plane, the magnification of each lens is 
constant, equal to 


M 0 = / d 0 (l) 

where d[ and d 0 are the image and object distances, respectively. Within the 
angular field indicated by the solid lines each lens images point objects with 
resolution and distortion that are within the specification of the lens. In the object 
plane the size of the region that can be imaged with good resolution by both 
lenses, the joint field’, is determined by the intersection of the respective angular 
fields. The lateral offsets d[^] and tf^are chosen to be the largest values that still 
permit the angular fields to overlap over the desired region in the object plane, i.e. 
to make the joint field as large as possible while still keeping the angle between 
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the lenses large enough to achieve small measurement error tor the ^-component ot 
velocity. Nominally, du and t/^2 are equal, and for the present system the angle 
between the axes ot the two lenses is 24.4 degrees. According to Prasad and Adrian 
( 1993) the error in the ^-component of velocity is about twice the error in the in- 
plane components of velocity when the angle between the axes of the lenses 
exceeds 30 degrees. 


JOINT FIELD 

OF VIEW OBJECT PLANE 



Figure 1 . Schematic of the stereoscopic camera. 
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The field of view of each CCD array is simply the area of the array 
mapped onto the object plane. As in Prasad and Adrian ( 1993) and Troy and Adrian 
( 1996), offsetting the cameras laterally in the image plane overlaps their respective 
fields of view. Each camera sees nearly the same region in the flow, thereby 
maximizing the size of the joint field of view. This design is relatively simple to 
lay out, but care must be taken to ensure that tolerances in the focal lengths of the 
lenses and in the orientations of the individual elements of the system are taken 
into account, either mechanically or by careful calibration. 

The joint field of view increases by more than a factor of two if the 
lenses are tilted towards each other so as to completely overlay their angular fields. 
This method, called ’angular displacement', has been discussed by Gauthier and 
Riethmuller (1988) and more recently by Prasad and Jensen (1996). Its principle 
advantage is the larger joint field of view, and its principal disadvantage is the 
more complicated alignment that is needed, which also involves using 
magnification that varies across the field of view. Only the simpler lateral offset 
method will be considered here, although angular offset can be easily 
accommodated by the present camera. 

The camera uses two Kodak Megaplus 42 CCD video cameras, each with 
2029 x 2044 pixels resolution. Since the camera is intended to image a 100-300 
mm field of view in the flow, and the CCD chip dimension is only 18 mm, the 
magnification of the system is substantially less than unity. M 0 = 0.23 for the 
present system. Two Nikon EL-Nikkor 302,5 mm F/5.6 enlarging lenses that 
are optimized for normal operation at large magnification are mounted in the 
reverse direction to provide the appropriate small magnification with little 
aberration, minimal distortion and large aperture (53.6 mm). The angular field of 
view of this lens is 57 degrees. Images are normally acquired at F/8. The video 
cameras and two lenses are mounted on a 3-D translation stage. The stage 
provides the six degrees of freedom needed to position and focus the cameras and to 
accurately register them with respect to each other. For the wind tunnel 
experiments, d Q = 1617 mm and d[ = 372mm. 

The rest of the stereoscopic system is shown in Fig, 2. Illumination of 
the flow field is provided by two Continuum Nd:YAG with a wavelength of 532 
nm at pulse frequency of 10 Hz. Each laser delivers up to 200 mJ/pulse of energy 
with a pulse duration of 8 nsec. The laser beam is formed into a sheet of about 1 
mm thickness in the test section. To resolve directional ambiguity, image 
shifting is normally needed. However image shifting was not used for the present 
measurements because of the strong mean flow in comparison with the turbulence 
intensity in the boundary layer. If directional ambiguity were a problem the 
cameras could be replaced by Ik x lk cross-correlating cameras. Timing for the 
sequence of laser pulses and the image capture with the cameras is provided by a 
TSI, Inc. Laser Pulse synchronizer box. 


3 Image Analysis 

The image analysis procedure is described in Figure 3. The capture and analysis of 
images is controlled by INSIGHT, a Windows-based software package supplied by 
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TSI, Inc. The software controls the simultaneous image capture by the video 
cameras. Image analysis to determine the velocity vectors on each image plane is 
performed by the autocorrelation technique. The outcome ot the analysis is a pair 
of displacement vector fields, AX] and AX 2 that are functions of the camera 
coordinates, X] and X 2 . The vectors in these fields are tested for validity using 
procedures described by Meinhart, et ai (1994), and invalid vectors are either 
replaced bv second or third choice vectors from the autocorrelation interrogation, 
replaced by interpolated values, or simply removed. Typically the Iraction ot 
invalid vectors is less than 2%. Corrections for the mis-registration ot the images, 
and the variation of the magnification factors are all performed by post-processing 
software. Effectively, this step maps the displacement vectors trom a grid in each 
of the CCD array image planes onto a common grid in the object plane. This step 
is necessary in general if there is any distorting medium between the object plane 
and the camera, such as a water-air interface (Prasad and Adrian 1993). In the 
present experiment the 6.35 mm thick glass windows of the wind tunnel cause the 
images to be displaced laterally by 0.55 mm. Since this shift is nearly constant, it 
is most easily accounted for by mechanical alignment in the initial layout of the 
cameras. In more general cases of nonlinear distortion it is better to calibrate the 
distortion and correct for it in software. The procedures for doing this will be 
discussed in another paper. 


doubt* tight ah*et 



Figure 2. Stereoscopic system. 
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Lastlw ihe stereoscopic equations are applied to compute ihe third- 
component ol the velocity vectors. For the ideal camera system shown in Fig. 1 
the image of a point particle at x is mapped by lens L f , onto a location rj f in the 

X p -plane on camera t given by 


-*J| = 



[(.r + <i tl )i + yy]. 


(2a) 


-Hz 



[{x-d L2 )i + yy} 


(2b) 


where carats denote unit vectors. Three-dimensional displacement of an image hv 
an amount Ax results in two-dimensional displacements Arj| and At)-, that are 

found by taking the increment of (2). The results depend on the c-location of the 
particle at the time of the first exposure. However, PIV correlation analysis gives 
an estimate of the volume average of the displacements of the particles that lie in 
the measurement volume defined by the thickness of the light sheet and the area of 
the interrogation spot. By integrating the equation for the image displacement over 
such a volume and by making use of the normally large ratio of the object distance 
to the light sheet thickness. Prasad and Adrian (1993) have shown that the volume 
averaged displacements (denoted by overbars) are given by 

AT), = M Q 


— . ( x+d L\) — 

Avh A: 

x + M 

— 1 — Az 

d o 

0 



Aq 2 


— (x-dr 9) — 

Ax+- ^-Az 

x + M 

A\+—Az 

\ 

0 

. ' d o _ 


(3b) 


Given the two two-dimensional vector displacements At^Arfe three ot 
the equations in (3) can be solved for Ax, Ay, Az, and the fourth can be used to 
obtain a second estimate of Ay . Alternatively, the two Av -equations can be added 
together to obtain a slightly more robust equation. The three-dimensional velocity 
vector is, of course, found by dividing three dimensional displacement by the time 
between exposures. 

Figures 4a and 4b show the displacement fields that result on cameras 1 
and 2, respectively when the object field undergoes a pure, uniform displacement 
in the ^-direction. Solving equation (3) for the three-dimensional particle 
displacement yields the field shown in Fig. 4c. 
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(a) 


(b) 



Figure 4. Image displacement fields due to a pure out-of-plane translation, (a) 
Camera 1, AXj(Xi); (b) AX 2 (X 2 ); (c) Three-dimensional displacement field 
Ax(X). 

4. Turbulent boundary layer measurements 

The turbulent boundary layer measurements were performed in an eiffel-type low- 
turbulence boundary layer wind tunnel with a working section 914 mm wide x 457 
mm high x 6.096 m long. The free stream turbulence intensity at the test section 
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inlet is less than 0.2% for free stream velocities less than 10 ms" 1 . The turbulent 
boundary layer was produced on a flat plate placed 100 mm above the floor of the 
test section. The boundary layer was tripped by the a 4.7 mm diameter wire which 
spanned the boundary layer plate just downstream of the leading edge. Optical 
access to the boundary layer was provided from the side by float glass windows, 
and from below by 610 mm wide x 2.748 m long float-glass windows embedded 
in the boundary layer plate. The stereo PIV measurements presented here were 
performed at a free stream velocity of 3.4 ms" 1 , which produced a Reynolds 
number based on the momentum thickness 9 of Re$= 2575 at the location of the 
measurements. All measurements were made with the light sheet in the x-y plane 
of the flow, where x is stream wise and y is normal to the wall. 

Figure 5 shows the (u,v) components of one realization of the flow field. 
A constant convection velocity of 3000 mm s* 1 has been subtracted to make the 


0.3 m/s 



Figure 5. Sample realization of the stream wise- wall normal f.v-vj vector field of 
the turbulent boundary layer. Re$— 2525. The free stream velocity is 3.4 ms' 1 , 

and a constant convection velocity of 3.0 ms" 1 has been subtracted to reveal the 
eddy structure more clearly. The vectors are the projection of the three-dimensional 
vectors onto the x-y plane. 
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fluctuations more visible. The edge of the turbulent boundary layer lies about 50 
mm above the wall in this realization. As observed in earlier work, there is an 
internal layer character close to the wall that grows in a manner similar to that of 
an ordinary boundary layer. This internal layer, seen as a region of uniformly low 
momentum that extends from the wail up to about 20 mm, is capped by a 
collection of several intense spanwise vortices. 

The w-component of velocity is shown in Figure 6 by plotting its 
contours in gray-level form. It is immediately clear that the values of the w- 
component are frequently as large or larger than the other components. There is a 
tendency for the w-contours to align along the same -30 degree inclination that is 
observed in the u-v field. There is also a tendency for the sign of the w-component 
to be positive on one side of a 30 degree line and negative on the other, indicating 
a large-scale rotation about an inclined axis. The strength of the w-component 
and its coherent organization shows that one must be very cautious in interpreting 
two-dimensional projections of the three-dimensional vectors such as the field in 
Figure 5. 



Figure 6. Sample realization of the streamwise-wall normal ( x-y ) vector field 
from Fig. 5 with the field of the spanwise turbulent velocity superimposed on it 
in the form of grey-level contours. The heavy line denotes the zero value of the 
w-component. and dashed lines denote negative values. 
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x(mm) 


Figure 7. Three-dimensional vectors on the .v-v plane in oblique 
view. All parameters as in Fig. 5. 



The entire plane of three-dimensional vectors is shown in oblique view in 
Figure 7. The growth of an internal low-momentum layer is very evident. The 
pattern is attributed to the alignment of inclined hairpin vortices (c.f. Meinhart and 
Adrian 1997, Zhou. Adrian and Balachandar 1997) In Figure 8 profiles of the \ -u 
components of the velocity vector are plotted for several different .v-locations. 
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These vectors demonstrate very clearly the oscillation in the sign of the w- 
component with increasing distance above the wall, and the inclination of the 
pattern at about thirty degrees with respect to the wall. The oscillation is 
consistent with the presence of inclined, hairpin-like vortices, as observed by 
Meinhart (1994) and Meinhart and Adrian (1997). 



Figure 8. Profiles of (v(x,y), w(x t y)) at several ^-locations from the vector field 
in Fig. 7. 

5 Summary 

Stereoscopic PIV corrects for perspective error as well as providing the out-of- 
plane component. It requires a special camera, but otherwise it uses essentially the 
same hardware as a monocular PIV. Videographic recording increases the ease of 
use substantially by eliminating the necessity of registering photographic images. 
The present system has been used to measure thousands of frames in experiments 
that required only a few days to set up. Although the video-based PIV system has 




The present system has been used to measure thousands of frames in experiments 
that required only a few days to set up. Although the video-based PIV system has 
less resolution than a photographic system, it is encouraging that it is still able to 
resolve many of the important features of the turbulent boundary layer that have 
only recently been discovered by photographic PFV. The results are very similar to 
those obtained by photographic PIV. Thus, while the video system clearly has a 
role to play in the acquisition of data for analysis by statistical averaging, it also 
provides enough resolution to be of considerable value in visualizing the flow. 


Acknowledgments 

This work was supported by NASA Grant NAG 1 1530, the United 
States Office of Naval Research, and TSI, Inc. 

References 


Adrian, R. J. 1991, Particle-imaging techniques for experimental fluid mechanics, 
Ann. Rev . Fluid Mech vol. 23, pp, 261-304. 

Arroyo, M.P. and Greater, C.A., 1991, Stereoscopic particle image velocimetry, 
Meas. Sci. Tech. Vol. 2, pp. 1181-1186. 

Barnhart, D.H., Adrian, R.J. and Papen, G.C., 1994, Phase-conjugate holographic 
system for high resolution particle image velocimetry, Appi Opt 
Vol. 33, pp. 7159-7170. 

Brodkey, R.S., 1977, Stereoscopic visual studies of complex turbulence shear 

flows, Flow Visualization, (Ed. Asanuma, T.), New York: McGraw-Hill, 
pp. 117-122. 

Brucker, C., 1995, 3D-PIV using stereoscopy and a scanning light-sheet: 

application to the 3D unsteady sphere wake flow. Flow Visualization 
VII , (Ed. Crowder, J.P.), New York: Begell House, pp. 715-720. 

Gauthier, V. and Riethmuller, M.R., 1988, Application of PIDV to complex 
flows: measurements of the third component. Von Karman Institute 
Lecture Series on Particle Image Velocimetry. Rhodes-Saint-Genese: 
von Karman Institute for Fluid Dymanics. 

Hinsch, K.D., 1995, Three-dimensional particle velocimetry, Meas y Sci. Tech. 
Vol. 6, pp. 741-753. 


14 



Meinhart, C.D.. 1994. Investigation of turbulent boundary-layer using particle 
image velocimetry. Ph. D. thesis. University of Illinois. 

Meinhart. C.D. and Adrian. R.J.. 1997, The coherent structure of the overlap 
region in a turbulennt boundary layer, (manuscript in preparation). 

Meng, H. and Hussain. F., 1995. In-line recording and off-axis viewing (IROV) 
technique for holographic particle velocimetry, Appl. Opt. Vol. 34. 
pp. 1827-1840. 

Nishino, Y„ Kasagi, S. and Hirata. M.. 1989, Three-dimensional particle tracking 
velocimetry based on automated digital image processing. ASME J. 

Vol. 111. pp. 384-391. 

Prasad. A.K.. and Adrian. R.J., 1993, Stereoscopic particle image velocimetry 
applied to liquid flows. Exp. Fluids. Vol. 15, pp. 49-60. 

Prasad, A.K., Jensen, K., 1995, Scheimpflug stereocamera for particle image 
velocimetry in liquid flows, Appl. Opt., Vol. 34, pp. 7092-7099. 

Robinson, O., and Rockwell, D., 1993, Construction of three-dimensional images 
of flow structure via particle tracking techniques, Exp. Fluids. Vol. 14, 
pp. 257-270. 

Troy, V. and Adrian R. J„ 1996, Investigation of turbulent penetrative 

convection, 1996, Eighth International Symposium on Applications of 
Laser Techniques to Fluid Mechanics , Lisbon, Portugal, 
pp.23.1.1-23.1.8.15 

Willert, C.E. and Gharib, M„ 1991, Digital particle image velocimetry, 

Exp. Fluids, Vol. 10, pp. 181-193. 

Zhou J., Adrian R. J„ Balachandar S. and Kendall T„ 1997, Hairpin vortices in 
near-wall turbulence and their regeneration mechanisms, (manuscript in 
preparation). 


15 



Report No. TP-95-101 


APPENDIX B 

Source Code for Combination of the Two-dimensional 
Displacement Fields from Stereoscopic PIV Images into Three- 

dimensional Vectors 

Z. C. Liu 



program stereo 


c **+***************+**+*+**************************************************** 

c This program is for combining the measurements of particle 

c displacements from two CCD cameras and to 

c calculate three components of velocity field. 

c The calculations are for the stereo-system with no window and no change 

c in the index of refraction in the optical paths, 

c The .p file assigns the basic parameters: 

c xregion, yregion: the mm size of the camera views 

c imax, jmax: the grid numbers in x and y directions 

c cmag: the average magnification of the two cameras 

c xshift, yshift: the px shifts of viewl (right) to viewO(left) 

c dt : the time separation of two exposures in seond 

c infileO, infilel: the .s files of cameraO and cameral, 

c from cleanvec output, x,y,u,v in px 

c outvec: output file of three velocity components in mm/s 

c interpolation: guassian interpolation 

£★★★+********★**+★★★★★★*★★★*★★+★**********★★+***★★*****★*★★****★*★***★**★**★★ 


include ' stereo. h' 

open (7 , f ile= ' stereo . p ' , status= ' unknown' ) 

read (7, ' (a) ' ) infileO 
read (7, ' (a) ' ) infilel 
read (7, ' (a) ' ) outvec 
read (7, ' (a) ' ) aorbin 
read (7,*) convect 
read ( 7 , * ) imax , j max 
read (7,*) xshif t , yshift 
read ( 7 , * ) cmag 
read ( 7 , * ) dt 

close (7) 

do=f * ( 1 . +1 . / cmag) 
di=f * (1 . +cmag) 
ms = cmag* s 

call data_in 

c xs and ys are in px, to convert them into grid numbers: 

xs=xshift/f loat (imax) 
ys=yshif t/ float (jmax) 

c To find the shifts of viewl to viewO in a grid unit: 

dx=xs-int (xs) 
dy=ys-int (ys) 

call gaussian 
call combine 

call vec out 


2000 stop 
end 



subroutine data in 


include 1 stereo.h' 
real xkern, ykern, choice 

c the input data are displacemants in px 

c multiply pxsize to convert px to mm in camera planes 

c the input data files are the .s files from cleanvec 

open (10, f ile=inf ileO, status= ’ unknown' , form= 1 unformatted' ) 

read (10) imax, j max, xkern, ykern 

read (10) choice , choice , choice , choice, choice 

do 200 j=l,jmax 
do 100 i=l,imax 

read (10) xO ( i , j ) , yO ( i , j ) , uO { i , j ) , vO ( i , j ) 
c write (*,*) xO ( i , j ) , yO ( i , j ) , uO ( i , j ) , vO ( i , j ) 

c The position data in .s file have already changed to mm, but the 

c velocity data are still pixels, not converted to mm not tnm/s 

xO (i, j ) =-x0 (i, j ) 
yO (i, j ) =-y0 (i , j ) 

c The bad vectors in .s file are given 9.e+6 values, and now changed 

c to zeros here 

c Also, if desire to look at flow field from a moving frame (ccd cameras) 

c with a convection velocity, subtract convect from u. The resulting three 

c components of velocity contain the difference only in u, while the v and 

c w do not change . 

if (uO ( i , j ) . eq . 9 . e+6 ) then 
uO ( i , j ) =0 . 0 
else 

uO (i , j ) =-u0 (i, j ) *pxsize 
endif 

if ( vO ( i , j ) . eq . 9 . e+6 ) then 
vO ( i , j ) =0 . 0 
else 

vO (i, j ) =-v0 (i, j ) *pxsize 
endif 

100 continue 

200 continue 

close (10) 

write (*,'(a)') infile0,'has been read.' 

open (15, file=infilel, status® 'unknown' , form® ' unformatted' ) 

read (15) imax, j max, xkern, ykern 

read (15) choice, choice, choice, choice, choice 

do 400 j=l,jmax 
do 300 i=l,imax 



read (15) xl (i, j ) ,yl (i, j ) ,ul (i, j ) 
c write (*,*) xl ( i , j ) , yl ( i , j ) , ul Ci , 

xl (i, j ) = - xl (i, j ) 
yl (i, j ) =-yl (i, j ) 

if (ul (i, j ) .eq. 9 .e+6) then 
ul (i , j ) =0 . 0 
else 

ul (i , j ) =-ul (i , j ) *pxsize 
endif 

if (vl (i, j ) .eq. 9 .e+6) then 
vl (i, j ) =0 . 0 
else 

vl (i, j ) =-vl (i, j ) *pxsize 
endif 


300 continue 

400 continue 

close (15) 

write (*, ' (a) ') infilel, 'has been read. ' 
write (*,*) 'data input is done!' 

return 

end 


, vl ( i , j ) 
j ) , vl (i, j ) 



subroutine gaussian 
include ' stereo. h' 


c the 2D-gaussian interpolation: 

nw * * ne 

I a | b | 

| (i, j ) | 

I I I 

c d 


sw * * S e 

c the values of x,y,u,v at (i,j) - the regular grid as the same as 

c the viewO - are the interpolated values at 4 points around (i, j) . 

c The interpolation is conducted by convolving the data with a gaussian 

c function. 

deltax=xO (2,1) -xO (1,1) 
deltay=yO (1,2) -yO (1,1) 

do 200 j=l,jmax 
do 100 i=l,imax 

inw=i+int (xs) 
ine=i+int (xs+1 . ) 
isw=inw 
ise=ine 

jnw= j + int (ys) 
jne= j +int (ys+l . ) 
j sw= jnw 
j se= jne 

a=dx*dx/deltax/deltax+dy*dy/ deltay/ deltay 

b= (deltax-dx) * (deltax-dx) /deltax/deltax+dy*dy/deltay/ deltay 
c=dx*dx/deltax/deltax+ (deltay-dy) * (deltay-dy) /deltay/ deltay 
d= (deltax-dx) * (deltax-dx) /del tax/de It ax+ 

+ (deltay-dy) * (deltay-dy) /deltay/ deltay 

a=exp ( -6 . *a) 
b=exp ( -6 . *b) 
c=exp ( -6 . *c) 
d=exp ( -6 . *d) 

xlp ( i , j ) = (xshif t+f loat (i*igrid) ) *pxsize 
yip (i , j ) = (yshift+f loat ( j* jgrid) ) *pxsize 

ulp (i, j ) =ul (inw, jnw) *a+ul (ine, jne) *b+ul (isw, jsw) *c+ul (ise, jse) *d 
vlp (i, j ) =vl (inw, jnw) *a+vl (ine, jne) *b+vl (isw, jsw) *c+vl (ise, jse) *d 

ulp ( i , j ) =ulp ( i , j ) / (a+b+c+d) 
vlp ( i , j ) =vlp ( i , j ) / (a+b+c+d) 

100 continue 

200 continue 


c 

c 

c 

c 

c 

c 

c 

c 

c 



write (*,*) 'Gaussian interpolation is done!' 

return 

end 



subroutine combine 


include 1 stereo. h' 

c convert registration shifts of cameral to cameraO into mm in camera 

c plane 

xsmm=xshi f t *pxs i ze 
ysmm=yshis t *pxs i ze 

c calculate the image region size into mm 

xregion= (xO (imax, 1) -xO (1,1) ) *pxsize 
yregion= (yO (1, jmax) -yO (1,1)) *pxsize 

c grid points (i,j) are in target plane, while the conjugate grid points c 

c the conversion of (xO,yO) in the image plane to (x,y) in the target 

c plane is being done in the following do loop 

c the combining formula have also taken care of the conversion of the c 

do 450 j=l,jmax 
do 440 i=l , imax 

ip=imax-i+l 
jp=jmax- j +1 

x ( i , j ) = (xregion-xO ( ip, jp) ) / cmag 
y (i, j ) = (yregion-yO (ip, jp) ) /cmag 

c if there are bad vector in (i,j) point of .s file, do not combine 

c and set them as zeros 

if (u0(ip,jp) . eq. 0 . 0 .or . vO (ip, jp) . eq . 0 . 0 . or . ulp ( ip, jp) .eq.0.0. 

+ or .vlp (ip, jp) .eq. 0 . 0) then 

w(i, j ) =0 . 0 
u ( i , j ) =0 . 0 
v ( i , j ) =0.0 

else 

w (i , j ) =do/ (1 . -ms/ (uO (ip, jp) -ul (ip, jp) ) ) 

u (i , j ) = (ulp (ip, jp) * (xO (ip, jp) -0.5* (ms+xregion) ) - 
+ uO (ip, jp) * (xlp (ip, jp) +0.5* (ms-xregion) -xsmm) ) / 

+ cmag/ (ms-uO (ip, jp) +ulp (ip, jp) ) 

v(i,j)=w(i,j)*0. 5/di* (vO (ip, jp) +vlp (ip, jp) - 
+ yregion+yO (ip, jp) +ylp (ip, jp-ysmm) ) - 

+ 0 . 5/cmag* (vO (ip, jp) +vlp (ip, jp) ) 

u ( i , j ) =u ( i , j ) /dt 
v ( i , j ) =v ( i , j ) /dt 
w(i,j)=w(i,j) /dt 

endif 

c write (*,*) x(i,j),y(i,j),u(i,j),v(i,j),w(i,j) 

440 continue 

450 continue 

write (*,*) 'Combination of the two views has been done!' 



return 

end 



subroutine vec out 


include ' stereo. h' 
b=0 . 0 

c c = c on ve c t *pxs i z e / d t / cmag 

c Two kinds of output formats can be selected: tecplot ascii with 

c a subtraction of a constant convection veclocity in streamwise 

c direction; or binary without any subtraction. 

300 if (aorbin . eq tecplot ' ) then 

open (30 , f ile=outvec , status= 1 unknown' , f orm= 1 formatted' ) 

write (30,*) 'TITLE = "BOUNDARY LAYER FLOW FIELD, Uo=3.4 m/s'" 
write (30,*) 'VARIABLES = X, Y, Z, U, V, W’ 

write (30,*) ' ZONE T= "CHANNEL_TEC" , 1= imax, ’ J= ', jmax, ' F=POINT ’ 

do 600 j=jmax,l,-l 
do 500 i=imax,l,-l 

if (abs (w { i , j ) ) . gt . 2000 . ) then 
w(i, j ) =0 . 0 
u(i, j ) =0 . 0 
v ( i , j ) =0 . 0 
else 
endif 

if (u (i, j ) .ne . 0 . 0) then 
u(i,j)=u(i,j) -cc 
else 
endif 

write (30,700) x ( i , j ) , y ( i , j ) , b, u ( i , j ) , v ( i , j ) , w ( i , j ) 

500 continue 

600 continue 

elseif ( aorb in. eq. ' binary ' ) then 

open (30 , f ile=outvec, status= ' unknown' , form= ' unformatted ' ) 

write (30) imax, jmax, xkern,y kern 

write (30) choice, choice, choice, choice, choice 

do 680 j=jmax,l,-l 
do 650 i=imax, 1,-1 


if (abs (w (i , j )) .gt . 2000 . ) then 
w ( i , j ) =0 . 0 
u (i, j ) =0 . 0 
v(i, j) =0.0 
else 
endif 

write (30) x ( i , j ) , y ( i , j ) , b , u ( i , j ) , v ( i , j ) , w ( i , j ) 


continue 

continue 


650 

680 


else 



+ 


700 


write (*,*) 'enter output f ormat ( tecplot or binary):' 
read (*,*) aorbin 

goto 300 
endif 

write (*,*) 'Calculation of three velocity components has been done 
write (*,*) 1 and stored in file: outvec 

write (*,*) 'The format of output is: ', aorbin 
write ( * , * ) 

write (*,*) 'The calculated velocity field is from a viewer with a 
constant convection velocity of : ' , cc, ' (mm/s) ' 
write ( * , * ) 

close (30) 

format (6fll.4) 

return 

end 



stereo.h 


parameter (m=128 , n=128 ) 
parameter (pxsize=0 . 009 ) 
parameter (s = 700 . 4 , f = 302 . 3 ) 

common xO , y 0 , uO , vO , xl , y 1 , ul , vl , xlp , yip , ulp , vlp , x , y , u , v , w , 
xshif t , yshif t , xregion, yregion, imax, jmax, dx, dy, dt , 
igrid, jgrid, cmag, cmagO , cmagl , do, di , ms , xs , ys , convect , 
inf ileO , inf ilel , outvec , interpolate, aorbin 

real xO (m, n) , yO (m, n) , uO (m, n) , vO (m, n) 
real xl(m,n) ,yl(m,n) ,ul(m,n) ,vl(m,n) 
real xlp(m,n) ,ylp(m,n) ,ulp(m,n) ,vlp(m,n) 
real x(m,n) ,y(m,n) ,u(m,n) ,v(m,n) ,w(m,n) 
real cmagO (m, n) , cmagl (m, n) , convect 

real do, di , ms , xshif t , yshif t , xs , ys , dx, dy, dt 
integer imax, jmax 

character*12 inf ileO , inf ilel , outvec, aorbin 
character *1 interpolate 



wt0_01 . s 
wtl_01 . s 
wt_01 
binary 
0.0 
40,80 
25 . , -9 . 

0.229 
0 . 00010 

stereo .p : 

input .s file from cleanvec for cameraO, binary 

input . s file from cleanvec for cameral, binary 

output file for 3 components, bin: wt_07, tecplot: wt_07.dat 

output format: tecplot (ascii) or binary 

convection velocity, Uc in px, 7.63 for wt_01, 11.45 for others. 

It is used only for tecplot output . 

This number was obtained from cleanvec 
imax,jmax -- grid number 

x and y registration shifts of cameral to cameraO, px 
averaged magnification of the stereo-camera system 
delta t, time separation of the two pulses, sec 
0.0001 for wt_0l 
0.00015 for others. 
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APPENDIX B 

Source Code for Combination of the Two-dimensional 
Displacement Fields from Stereoscopic PIV Images into Three- 

dimensional Vectors 

Z. C. Liu 



program stereo 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 




**********************************************************+*+********* 

This program is for combining the measurements of particle 
displacements from two CCD cameras and to 
calculate three components of velocity field. 

The calculations are for the stereo-system with no window and no change 
in the index of refraction in the optical paths . 

The .p file assigns the basic parameters: 

xregion, yregion: the mm size of the camera views 
imax, jmax: the grid numbers in x and y directions 
cmag: the average magnification of the two cameras 
xshift, yshift: the px shifts of viewl (right) to viewO(left) 
dt : the time separation of two exposures in seond 
infileO, infilel: the .s files of cameraO and cameral, 

from cleanvec output, x,y,u,v in px 

outvec: output file of three velocity components in mm/s 
interpolation: guassian interpolation 




include 1 stereo . h ' 

open (7 , file® 1 stereo. p 1 , status® ' unknown' ) 

read (7, ' (a) ' ) infileO 
read (7, ' (a) ' ) infilel 
read (7, ' (a) ') outvec 
read (7, ' (a) ’) aorbin 
read ( 7 , * ) convect 
read ( 7 , * ) imax , j max 
read (7,*) xshif t , yshift 
read ( 7 , * ) cmag 
read ( 7 , * ) dt 

close (7) 

do=f * (1 . +l . / cmag) 
di = f * (l . -t-cmag) 
ms=cmag*s 


call data_in 

xs and ys are in px , to convert them into grid numbers : 

xs®xshift/float (imax) 
ys=yshift/float (jmax) 

To find the shifts of viewl to viewO in a grid unit: 

dx=xs-int (xs) 
dy=ys-int (ys) 

call gauss ian 
call combine 
call vec out 


2000 stop 
end 



subroutine data in 


include ' stereo. h' 
real xkern, ykern, choice 

c the input data are displacements in px 

c multiply pxsize to convert px to mm in camera planes 

c the input data files are the .s files from cleanvec 

open (10, f ile=inf ileO , status= 'unknown' , form= 'unformatted' ) 

read (10) imax, j max, xkern, ykern 

read ( 10 ) choice, choice, choice, choice, choice 

do 200 j=l,j max 
do 100 i=l,imax 

read (10) x0 (i , j ) , y0 ( i , j ) , uO (i, j ) , vO (i , j ) 
c write (*,*) xO (i , j ) , yo ( i , j ) , uO (i, j ) , vO (i , j ) 

c The position data in .s file have already changed to mm, but the 

c velocity data are still pixels, not converted to mm not mm/s 

x0 (i, j ) =-x0 (i, j ) 
yO (i, j ) =-y0 (i, j ) 

c The bad vectors in .s file are given 9. e+6 values, and now changed 

c to zeros here 

c Also, if desire to look at flow field from a moving frame (ccd cameras) 

c with a convection velocity, subtract convect from u. The resulting three 

c components of velocity contain the difference only in u, while the v and 

c w do not change . 

if (uO ( i , j ) . eq . 9 . e+6 ) then 
uO (i, j ) =0 . 0 

else 

uO (i, j ) =-u0 (i , j ) *pxsize 
endif 

if ( vO ( i , j ) . eq . 9 . e+6 ) then 
vO (i, j ) =0 . 0 

else 

vO (i, j ) =-v0 (i, j ) *pxsize 
endif 

100 continue 

200 continue 

close (10) 

write (*, ' (a) ' ) infileO, 'has been read. ’ 

open ( 15 , f ile=inf ilel , status= ' unknown ' , form= ' unformatted ' ) 

read (15) imax, j max, xkern, ykern 

read (15) choice, choice, choice, choice, choice 

do 400 j =1 , j max 
do 300 i=l , imax 



c 


read (15) xl (i, j ) , yi (i, j ) /Ul (i> j ) , vl (i , j ) 
write (*,*) xl ( i , j ) , yl ( i , j ) , ul Ci , j ) , vl ( i , j ) 


xl ( i , j ) =-xl ( i , j ) 
yl (i , j ) =-yl (i, j ) 

if (ul ( i , j ) . eq . 9 . e+6 ) then 

ul (i, j ) = 0.0 

else 

ul (i , j ) =-ul (i, j ) *pxsize 
end if 


if (vl ( i , j ) . eq . 9 . e+6 ) then 

vl (i, j ) =o . 0 

else 

vl (i, j ) =-vl (i , j ) *pxsize 
endif 


300 continue 

400 continue 

close (15) 

write (*, ' (a) ' ) infilel, 'has been read. 1 
wr i^ e (*,*) 'data input is done!' 

return 

end 



subroutine gaussian 

include ' stereo. h' 


the 2D-gaussian interpolation: 
nw * * ne 

I a | b I 

! 1 i 

j (i, j ) j 

! 1 1 

I C d I 


sw * * se 

the values of x,y,u,v at (i,j) - the regular grid as the same as 
the viewO - are the interpolated values at 4 points around (i,j). 

The interpolation is conducted by convolving the data with a gaussian 
function. 

deltax»xO (2,1) -xO (1,1) 
deltay-yO (1,2) -yO (1, 1) 

do 200 j = 1 , j max 
do 100 i=l,imax 

inw=i+int (xs) 
ine=i+int (xs+l . ) 
isw=inw 
ise=ine 

jnw= j+int (ys) 
jne= j+int (ys+l . ) 
j sw= jnw 
j se= jne 

a=dx*dx/ deltax/deltax+dy*dy/ deltay/deltay 
b= (deltax-dx) * (deltax-dx) /deltax/deltax+dy*dy/deltay/deltay 
c=dx*dx/deltax/deltax+ (deltay-dy) * (deltay-dy) /deltay/deltay 
d= (deltax-dx) * (deltax-dx) /deltax/deltax+ 

(deltay-dy) * (deltay-dy) /deltay/deltay 

a-exp(-6.*a) 
b«exp ( -6 . *b) 
c=exp ( -6 . *c) 
d=exp ( -6 . *d) 

xip ( i , j ) = (xshift+f loat ( i*igrid) ) *pxsize 
yip (i , j ) = (yshift + f loat ( j *jgrid) ) *pxsize 

ulp ( i , j ) =ul ( inw, jnw) *a+ul (ine, jne) *b+ul (isw, jsw) *c+ul (ise, jse) *d 
vlp (i, j ) =vl (inw, jnw) *a+vl (ine, jne) *b+vl (isw, jsw) *c+vl (ise, jse) *d 

ulp ( i , j ) =ulp ( i , j ) / (a+b+c+d) 
vlp ( i , j ) = v lp ( i , j ) / (a+b+c+-d) 


continue 

continue 



wr ^ e (*/*) 'Gaussian interpolation is done!’ 

return 

end 



subroutine combine 


include ' stereo. h' 

convert registration shifts of cameral to cameraO into mm in camera 
plane 

xsmm=xshi f t *pxs i ze 
ysmm=yshis t *pxs i ze 

calculate the image region size into mm 

xregion* (xO (imax, 1) -xO (1,1)) *pxsize 
yregion- (yO (1, jmax) -yO (1,1) ) *pxsize 

grid points (i,j) are in target plane, while the conjugate grid points 
the conversion of (xO,yO) in the image plane to (x,y) in the target 
plane is being done in the following do loop 

the combining formula have also taken care of the conversion of the 

do 450 j =1 , jmax 
do 440 i=l , imax 


ip=imax-i+l 
jp=jmax- j+1 

x (i , j ) = (xregion-xO (ip, jp) ) /cmag 
y (i, j ) = (yregion-yO (ip, jp) ) /cmag 

if there are bad vector in (i,j) point of .s file, do not combine 
and set them as zeros 

if (u0(ip, jp) . eq. 0 . 0 .or . v0 (ip, jp) . eq. 0 . 0 . or . ulp (ip, jp) . eq.0.0. 

or . vlp ( ip , jp) .eq.0.0) then 
w(i, j ) =o .0 
u(i, j)=0.0 
v(i, j) -o.o 

else 

w(i , j ) =do/ (1 . -ms/ (uO (ip, jp) -ul (ip, jp) ) ) 

u (i , j ) = (ulp (ip, jp) * (xO (ip, jp) - o . 5* (ms+xregion) ) - 

uO (ip, jp) * (xlp (ip, jp) +o . 5* (ms-xregion) -xsmm) ) / 
cmag/ (ms-uO (ip, jp) +ulp (ip, jp) ) 

v(i,j)=w(i,j) *0 . 5/di*(v0 (ip, jp) +vlp (ip, jp) - 
yregion+yO (ip, jp) +ylp (ip, jp-ysmm) ) - 
0 . 5/cmag* (vO (ip, jp) +vlp (ip, jp) ) 

u(i,j)=u(i,j) /dt 
v(i, j) =v(i, j ) /dt 
w(i,j)=w(i,j) /dt 

endif 

write (*,*) x(i,j) ,y(i,j) ,u(i,j) ,v(i,j) , w ( i , j ) 

continue 

continue 

write (*,*) 'Combination of the two views has been done!' 



return 

end 



c 

c 

c 

300 


500 

600 


650 

680 


subroutine vec_out 

include ' stereo. h' 
b=0 . 0 

cc=convect*pxsize/dt/cmag 

Two kinds of output formats can be selected: tecplot ascii with 
a subtraction of a constant convection veclocity in streamwise 
direction; or binary without any subtraction. 

if (aorbin. eq. 1 tecplot ' ) then 

open (30 , f ile=outvec , status= ' unknown ' , f orm= 1 formatted ' ) 

write (30,*) 'TITLE = "BOUNDARY LAYER FLOW FIELD, Uo=3.4 m/s" ' 
write (30,*) 'VARIABLES = X, Y, Z, U, V, W' 

write (30,*) ' ZONE T= "CHANNEL_TEC" , 1= 1 , imax, ' J* ' , jmax, ' F=POINT' 

do 600 j=jmax,l,-l 
do 500 i=imax, 1,-1 

if (abs (w(i, j ) ) .gt.2000.) then 
w(i , j ) =0 . 0 
u(i , j ) =0 . 0 
v ( i , j ) =o . 0 
else 
end if 

if (u(i, j) .ne.0.0) then 
u(i,j)=u(i,j) -cc 
else 
end if 

write (30,700) x ( i , j ) , y ( i , j ) , b, u ( i , j ) , v ( i , j ) , w ( i , j ) 


continue 

continue 

elseif (aorbin. eq. 1 binary' ) then 

open (30 , f ile=outvec , status* ' unknown ' , form* ' unformatted ' ) 

write (30) imax, jmax, xkem,ykem 

write (30) choice, choice, choice, choice, choice 


do 680 j*jmax, 1, -1 
do 650 i*imax,l,-l 

if (abs(w(i,j) ) .gt.2000. ) then 
w(i, j ) =0 . 0 
u (i , j ) =0 . 0 
v ( i , j ) =0 . 0 
else 
endif 

write (30) x (i , j ) , y ( i , j ) , b, u ( i , j ) , v (i , j ) , w ( i , j ) 


continue 
cont inue 

else 



write {*,*) 'enter output format ( tecplot or binary) 
read (*,*) aorbin 
goto 300 
endif 


+ 


write ( * , *) 
write (*,*) 
write (*,*) 
write (*,*) 
write (*,*) 
constant 
write (*, *) 


'Calculation of three velocity components has been done 
' and stored in file: ' , outvec 

'The format of output is: ' , aorbin 

'The calculated velocity field is from a viewer with a 
convection velocity of : 1 , cc, ' (mm/s) ' 


close (30) 


700 format (6fll.4) 
return 
end 



stereo . h 


parameter (m=128 , n=128 ) 
parameter (pxsize=0 . 009) 
parameter (s=700 . 4 , f = 302 . 3 ) 

common xO , yO , uO , vO , xl , yl , ul , vl , xlp, yip, ulp, vlp, x, y , u, v, w, 
xshif t , yshif t , xregion, yregion, imax, jmax, dx, dy, dt , 
igrid, jgrid, cmag, cmagO , cmagl , do, di , ms, xs , ys , convect , 
inf ileO , inf i lei , outvec , interpolate , aorbin 

real xO (m,n) ,y0 (m,n) ,u0 (m,n) , vO (m,n) 
real xl(m,n) ,yl(m,n) ,ul(m,n) ,vl(m,n) 
real xlp (m, n) , yip (m, n) , ulp (m, n) , vlp (m, n) 
real x(m,n) ,y(m,n) ,u(m,n) ,v(m,n) ,w(m,n) 
real cmagO (m, n) , cmagl (m, n) , convect 

real do, di , ms , xshif t , yshif t , xs , ys , dx, dy, dt 
integer imax, jmax 

character* 12 inf ileO , inf i lei , outvec, aorbin 
character *1 interpolate 



wt0_01 . s 
wtl_01 . s 
wt_01 
binary 
0 . 0 
40,80 
25. , -9. 

0.229 
0 . 00010 

stereo. p: 

input .s file from cleanvec for cameraO, binary 

input .s file from cleanvec for cameral, binary 

output file for 3 components, bin: wt_07, tecplot: wt_07.dat 

output format: tecplot (ascii) or binary 

convection velocity, Uc in px, 7.63 for wt_01, 11.45 for others. 
It is used only for tecplot output. 

This number was obtained from cleanvec 
imax,jmax -- grid number 

x and y registration shifts of cameral to cameraO, px 
averaged magnification of the stereo-camera system 
delta t, time separation of the two pulses, sec 
0.0001 for wt_01 
0.00015 for others. 



